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Abstract. In the Earth auroral zone, the electron acceleration by Alfven waves is sometimes a precursor 
of the non-propagating acceleration structures. In order to investigate how Alfven waves could generate 
non-propagating electric fields, a series of simulations of counter-propagating waves in a uniform plasma is 
presented. The waves (initially not configured to accelerate particles) propagate along the ambient magnetic 
field direction. It is shown that non propagating electric fields are generated at the locus of the Alfven waves 
crossing. These electric fields have a component orientated along the direction of the ambient magnetic 
field, and they generate acceleration and a significant perturbation of the plasma density. The non-linear 
interaction of down and up-going Alfven waves might be a cause of plasma density fluctuations (with gradients 
along the magnetic field) on a scale comparable to those of the Alfven wavelengths. 
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1 Introduction 



An auroral substorm is an abrupt increase in night-side auroral power. The largest part of the electron accel- 
eration that triggers the auroras happens at a few thousands of kilometres above the ionosphere. Two main 
families of acceleration processes have been identified: those associated to quasi-static electric fields called strong 
double layers, and those associated to Alfven wave electric fields (see ( |Mottez||2012b| ) for a review on the role of 
the Alfven waves). The acceleration by quasistatic, and basically non-propagating, electric fields such as double 
layers produces mono-energetic beams of electrons, while those by Alfven waves are associated to broadband 
energy distributions. Recent observations show that these two families of processes are not independent from 
each other. 

Evidences of acceleration structures emanating from Alfven waves are given by the direct observation of the 



parallel electric field (parallel to the mean direction of the magnetic field) (Chust et al. 



2007D or by the estimate of the wave Poynting flux ( |Louarn et al.|1994|[Volwerk et al.|1996)|Keiling et al.|2QQQ^ 



1998 Chaston et al. 



Several papers suggest that Alfvenic processes might act as the precursors of quasi-static non-propagating 
acceleration structures QZou et al.||2010| |Newell et al.||2010[ |Hull et al.||2010 ). 

This paper summaries a series of numerical simulations that investigate the interaction of down-going incident 
Alfven waves with up-going Alfven waves reflected on the ionosphere. Their ability to create parallel stationary 
electric fields is questioned, as well as to prepare the auroral plasma before the setting of the acceleration 
processes. 

In ideal MHD, Alfven waves do not carry a parallel electric field. Therefore, they cannot accelerate auroral 
electrons. It has been shown in previous studies (Hase gawa fc Chen||1975 Goertz||1984 ) that Alfven waves with 
an oblique wave vector induce an electric field with a component parallel to the ambient magnetic field. In the 
present study, we don't suppose that the Alfven waves are already accelerating the plasma. As the origin of 
their small transverse scales is not trivial, it is more careful, for the generality and the simplicity of the initial 
condition to neglect the presence of small transverse scale associated to the Alvfen waves. 

Because of the simple setting of the simulations, the present study cannot pretend to bring any conclusion on 
the large scale structure of the auroral zone. For instance, we cannot involve the large scale density and magnetic 
field dependence with altitude, net potential drops on large scales, etc. This paper focuses on an explanation 
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of the physical process observed in the simulations that are in the range of "microphysics", when compared to 
magnet ospheric scales. Nevertheless, its potential relevance to auroral physics is the main motivation of this 
work. 

Let us notice already that we make a distinction between the X component of a vector (parallel to the 
ambient -and uniform- magnetic field) and its x component, that is parallel to the local value of the magnetic 
field. (The reason why is explained in section |4|) 



2 Numerical method and simulation parameters 



As in flGenot et aT]|2000| |Genot et aL]|2001| |Genot et aL]|2004| |Mottez fe Genot||201l| ), dedicated to the 



physics of auroral Alfven waves, the numerical simulations are made with a EGC (Electron Guiding Center) 
electromagnetic PIC code that takes into account the motion of the electron guiding center, and the full ion 
motion. The boundary conditions are periodic. A complete description of this code is given in ( |Mottez et al.| 
1998). The method of initalisation of the Alfven waves is provided in (Mottez 2008), and the method for the 



wave packets in |Mottez| ( |2012a| ) 

The physical variables are reduced to dimensionless variables. Time and frequencies are normalized by 
the electron plasma frequency uj p q that correspond to a reference background electron density no- Velocities 
are normalized to the speed of light c, and the magnetic field is given in terms of the dimensionless electron 
gyrofrequency uj ce /uj p Q. The mass unit is the electron mass m e . Therefore, the units (starting from the Maxwell 
Eq. in the MKSA system) are c/uj p0 for distances, ou p o/c for wave vectors, e for charges, en for the charge 
density, cuj ce /uu p o for the electric field. In the following parts of this paper, all the numerical values and figures 
are expressed in this system of units. 



3 Simulations 



The left-hand side of Fig.l shows the magnetic field Bz of two wave packets that propagate in opposite direction. 
The field Bz is a good proxy of the wave packets positions. We can see that they cross each other at x ~ 120 and 
their intersection starts at t = 200. The right-hand side of Fig 1 shows the parallel electric Ex field associated 
to the same wave packets. The alternating fine horizontal lines are associated to the plasma oscillations that 
are present in any plasmas (their frequency is oj p o). Apart from the plasma waves, we can see that, as predicted 
by MHD laws, Ex is (almost) null before the intersection. But during the wave packets crossing, it is strongly 
enhanced. After the intersection, with a small delay, Ex becomes null again. In order to investigate this 
phenomenons, the problem has been simplified agin. Instead of two wave packets, the crossing of two sinusoidal 
waves has been studied. The resulting parallel electric field Ex is shown on Fig. 2. Here again, we can see the 
plasma waves, but also a time independent structure. Because it is present at the start of the simulation, it is not 
the consequence of an instability (it would not grow instantaneously). Actually, with the two monochromatic 
waves, the waves interaction starts at the beginning of the simulation since the two waves are present everywhere 
with the same intensity. 



4 How the parallel electric field sets 



An analysis of the properties of Ex was conducted in Mottez (2012a), and an explanation was provided. First, 



it was shown that the intensity of the field is proportional to the product of the intensity of the two sinusoidal 
waves. Then, it was shown that Ex is always independent of time but dependent on X, or independent of X 
and dependent on time. Then, it was shown that there is no critical value for this phenomenon to appear, it 
occurs also with low intensity waves, until it is below the noise level. We also projected the electric field along 
the local and instantaneous direction x of the magnetic field (that is not uniform because of the perturbation 
5Bi and 5B2 of the two waves). We found that E x — 0. Then, it was shown that it could be explained in the 
following way. The electric field is assumed to be perpendicular to the direction of the local magnetic field: 
E • B = 0. This property is true in the linear theory, it is generalized to non- linear interactions ( |Knud sen 1996; 
Tsiklauri 2007). Therefore, E x — 0, and the projection of the resulting electric field on the X axis can be 



derived. This derivation (provided in |Mottez| ( |2012b| )) depends on the polarization of each wave (righ-handed 
of left-handed circular polarization). It appears that Ex depends of products of EiE 2 , of E\ and of E\ where 
Ei and E 2 are the amplitudes (set initially) of the two waves. A second hypothesis is made at this time: we 
consider that the rule E • B = is true as long as it concerns the interaction of two different waves, but not 
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Fig. 1. Plots of electromagnetic field components as a function of position X (horizontal axis) and time t (vertical 
axis).Le/i- The magnetic field Bz associated to two wave packets that propagate in opposite direction. Right: the 
parallel electric field Ex associated to the same wave packets. 



F_p Xpar Vpar Time 0.00 




-0.130102E-01 0.129180E-01 0.00 22.00 



Fig. 2. Left: The parallel electric field Ex(X,t) associated to a pair of monochromatic Alfven waves that propagate in 
opposite directions. (Same axis as for Fig. 1). Right: The electron phase space density for the same simulation as for 
Fig. 1. The horizontal axis represents the direction X (it is parallel to the ambient magnetic field), the vertical axis 
represents the parallel electron velocity Vx and the grey scale is the electron density in that space. 
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of a wave with itself. (Why ? This point still needs to be clarified.) Then, in the expression of Ex, only the 
terms proportional to EiE 2 are kept. They perfectly match the properties found in the simulations analysis, 
according to the various choices of amplitudes, directions of propagation and polarisations. 

5 Particle acceleration and relevance to auroral physics 

The previous analysis shows that the various waves that contribute to a wave packet also interact and they 
contribute to a parallel electric field (even for a single wave packet). This is why Ex is not strictly null at the 
begining of the simulation. The figure 3 shows the phase space of the electrons at three different times in the 
same simulation as for Fig 1. We can see that there is initially a perturbation of the parallel velocity of the 
electron associated to the initial parallel field Ex- 

More interestingly, we can see that after the wave crossing, the parallel velocity is locally strongly enhanced, 
and this enhancement lasts well after the wave packets crossing, contrarily to the field Ex- This electron 
distribution in the phase space presents interesting similarities with those of a newly settle strong double layer 
(localized strong acceleration, well above the thermal level, with a shift of the bulk electron distribution without 
heating, thus ready to provide quasi mono-energetic electron beams). This is interesting, because as it was said 
in the introduction, these electrostatic structures (strong double layers) dominate the auroral acceleration after 
the phase of the Alfvenic processes. 

6 Conclusion 

Two Alfven wave packets crossing each other generate an electric field in a direction X that is parallel to the 
average ambient magnetic field. This can be explained if we consider that the two waves interact in a way 
that let the wave electric field and the total magnetic field perpendicular to each other. This is a non-linear 
wave- wave interaction whose intensity is characterized by the intensity of each wave. The electric field E x is 
favourable to electron acceleration, and the phase-space distribution of the electron keeps a signature of the two 
waves interaction well after the waves crossing has occurred. The influence of the accelerated electron will be 
the object of a further study. 

The simulations presented here may provide important clues explaining the transition from the Alfvenic to 
the electrostatic auroral acceleration processes mentioned in the introduction. 
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